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ABSTRACT 

We present a study of extended galaxy halo gas through HI and OVI absorption 
over two decades in projected distance at z « 0.2. The study is based on a sample 
of 95 galaxies from a highly complete (> 80%) survey of faint galaxies (L > O.IL*) 
with archival quasar absorption spectra and 53 galaxies from the literature. A clear 
anti-correlation is found between HI (O VI) column density and virial radius nor¬ 
malized projected distance, d/Rh- Strong HI (OVI) absorption systems with column 
densities greater than 10^^ ° (10^^'®) cm“^ are found for 48 of 54 (36 of 42) galaxies 
at d < i?h indicating a mean covering fraction of (khi) = 0.89 ((kqvi) = 0.86). 
OVI absorbers are found at d « i?h, beyond the extent observed for lower ionization 
species. At d/Rh = 1 i strong HI (OVI) absorption systems are found for only 7 
of 43 (5 of 34) galaxies ((khi) = 0.16 and (kqvi) = 0.15). Beyond d = 3i?h, the HI 
and OVI covering fractions decrease to levels consistent with coincidental systems. 
The high completeness of the galaxy survey enables an investigation of environmen¬ 
tal dependence of extended gas properties. Galaxies with nearby neighbors exhibit a 
modest increase in OVI covering fraction at d > i?h compared to isolated galaxies 
(kqvi ~ 0-13 versus 0.04) but no excess HI absorption. These findings suggest that 
environmental effects play a role in distributing heavy elements beyond the enriched 
gaseous halos of individual galaxies. Finally, we find that differential HI and OVI 
absorption between early- and late-type galaxies continues from d < i?h to d « 3i?h- 

Key words: galaxies: haloes - quasars: absorption lines - galaxies: interactions 


1 INTRODUCTION 

Galaxies grow by a combination of gas accretion from 
the intergalactic medium (IGM) and mergers. After heavy 
element enrichment, it is expected that gas is subsequently 
ejected back to the IGM through starburst and AGN driven 


winds (e.g. Murray et al. 20111 or by stripping during 


galaxy interactions (e.g. Wang 1993 Gauthier 20131. The 
low-density gas of the circum-galactic medium (GGM) is at 
the nexus of these baryon cycles, and it represents a po¬ 
tentially significant, multiphase reservoir that can fuel star- 
formation on long time-scales (e.g. Mailer &: Bullock|2004 1. 
Moreover, the CGM constitutes a sensitive laboratory for 
studying the physical processes that control the inflow and 
outflow of baryons from galaxies and for testing sub-grid 
feedback recipes implemented in simulations of galaxy evo¬ 
lution (e.g. Hummels et al.|2013 Ford et al.|2013 Gen|2013 


Shen et al. 2013 Agertz fc Kravtsov||2014). 

The density of the GGM is nearly always too low to be 


studied in emission with existing facilities. Nevertheless, a 
great deal of progress has been made via UV absorption¬ 
line spectroscopy of background objects at low projected 
distances, d, from foreground galaxies. Observations of the 


HI Lyman series (e.g. 

Chen et al. 

1998 [Tripp et al. 

1998 

Wakker & Savage|2009 

[Stocke et a 

.|2013||Rudie et al. 

2013 


the MgII doublet (e.gJBowen(et^^ 

1995 |Ghen et al'|201^ [Gauthier et al.|2010| [Bordoloi et al. 


,2011 1, the GIV doublet (e.g. Chen et al. 2001 Borthakur 
et al. 2013 Liang & Chen 2014 Bordoloi et al. 20141, and the 


O VI doublet (e.g. Chen & Mulchaey 2009 Wakker & Savage 
2009[ Prochaska et al.||2011 Tumlinson et al.||2011||Mathes 


et al. 2014 Stocke et al. 2014 ?) have been particularly 


fruitful. 

Cool, metal-enriched gas traced by MgII, CII, Sill, 
Si III, and CIV absorption is observed out to impact pa¬ 
rameters of 0.7 times the galaxy virial radius (d ~ 0.7i?h) 
around both late-type and early-type galaxies of L < L* 
alik(0 At larger projected distances, the incidence of these 


E-mail: seanjohnson@uchicago.edu 


^ There is evidence that CIV absorption is more extended around 
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Table 1. Summary of galaxy and absorber properties. 


Quasar 

sightline 

Galaxy 

R.A. Dec. 

(J2000) (J2000) 

2gai log (M^/M0)c,r class^ 

envi.'^ 

survey 

d 

(kpc) 

d/Rh 

log A/cm ^ 

HD 0 VD 

J0042-1037 

00:42:22.27 

-10:37:35.2 

0.0950 

9.5 

Late 

I 

COS-Halos 

15 

0.1 

14.80 - 18.50 

14.70 + 0.22 

J0226±0015 

02:26:12.98 

+00:15:29.1 

0.2274 

10.4 

Early 

I 

COS-Halos 

81 

0.4 

14.36 + 0.06 

< 13.12 

HE0226-4110 

02:27:46.27 

-40:53:15.1 

0.1120 

10.4 

Early 

I 

IMACS 

826 

4.1 

< 12.88 

< 12.90 

HE0226-4110 

02:27:57.47 

-40:57:25.7 

0.2952 

8.6 

Late 

I 

IMACS 

884 

9.8 

13.36 + 0.05 

< 12.75 

HE0226-4110 

02:28:01.44 

-40:57:55.7 

0.3863 

9.8 

Late 

I 

IMACS 

846 

6.5 

13.39 + 0.04 

< 12.66 

HE0226-4110 

02:28:04.31 

-40:55:57.4 

0.1772 

10.5 

Early 

NI 

IMACS 

435 

2.1 

< 12.23 

< 13.07 

HE0226-4110 

02:28:04.52 

-40:57:02.3 

0.3971 

11.1 

Late 

NI 

IMACS 

648 

1.8 

13.99 + 0.02 

13.57 + 0.09 

HE0226-4110 

02:28:09.37 

-40:57:44.5 

0.3341 

9.7 

Late 

NI 

IMACS 

346 

2.6 

< 12.36 

< 12.75 

HE0226-4110 

02:28:10.23 

-40:55:48.1 

0.1248 

8.5 

Late 

I 

IMACS 

230 

2.3 

12.88 + 0.05 

< 12.84 

HE0226-4110 

02:28:10.50 

-40:58:57.3 

0.3268 

9.4 

Early 

I 

IMACS 

547 

4.7 

< 12.35 

< 13.05 

HE0226-4110 

02:28:11.32 

-40:56:36.7 

0.2804 

8.4 

Late 

I 

IMACS 

244 

2.8 

13.73 + 0.01 

-1 

HE0226-4110 

02:28:11.59 

-40:58:53.8 

0.2492 

8.6 

Late 

I 

IMACS 

419 

4.5 

13.38 + 0.03 

< 13.16 

HE0226-4110 

02:28:13.89 

-40:59:46.2 

0.3939 

10.2 

Late 

I 

IMACS 

812 

5.2 

< 12.36 

< 12.62 

HE0226-4110 

02:28:14.56 

-40:57:22.7 

0.2065 

8.9 

Late 

NI 

IMACS 

37 

0.4 

15.27 + 0.06 

14.37 + 0.01 

HE0226-4110 

02:28:16.29 

-40:57:27.2 

0.2678 

10.3 

Early 

I 

IMACS 

75 

0.4 

< 12.97 

< 12.87 


Notes 

The full table is available in the online version of the paper. 

^ Galaxy classification with “Early” for early-type, absorption-line dominated galaxies and “Late” for late-type, emission-line dominated galaxies 
^ Galaxy environment class as defined in the text with “I” for isolated, “NI” for non-isolated, and “A” for ambiguous cases. 

^ Total column density measured in the COS sightline within Av = ±300 kms~^ of the galaxy systemic redshift. In cases where O VI or HI column 
densities cannot be measured due to contamination, the value is set to -1. 


absorption species decreases sharply to near zero. |Liang| 
& Chen (20141 present two possible explanations for this 


“metal-boundary”: ( 1 ) an inability of galactic outflows and 
satellite accretion to efficiently deposit heavy element en¬ 
riched gas at distances greater than 0.7 Rh or (2) an inabil¬ 
ity of cool-warm gas clouds to form and survive at larger 
distances. 

In contrast, highly-ionized and enriched gas traced by 
the OVI doublet is found to reflect the bimodal distribu¬ 
tion of galaxies in color-magnitude diagrams. In particu¬ 
lar, strong OVI systems are found at high incidence out to 
d = 150 kpc around star-forming galaxies but at lower inci¬ 
dence around more quiescent ones ( |Chen &: Mulchaey|2009[ 
Tumlinson et al. 2011 1 . Moreover, OVI detections have been 


reported around sub-L* galaxies out to d = 300 kpc which 
corresponds to 1 — 3i?h ( |Prochaska et alT]|201ip . However, 
the sample of galaxies probed at such large distances by ab¬ 
sorption spectroscopy is small and the spatial extent of O VI 
bearing gas around galaxies remains an open question. 

To examine the spatial dependence of OVI absorption 
around galaxies, we have assembled a sample of 95 galax¬ 
ies that are probed in absorption by quasar sightlines with 
UV spectra from the Cosmic Origins Spectrograph (COS; 
Green et al. 20121 on the Hubble Space Telescope (HST). 
These galaxies are drawn from a highly complete (> 80% 
complete for galaxies of L > O.IL* at z < 0.4 and d < 500 
kpc) survey of faint galaxies in the fields of four high signal- 
to-noise COS quasar sightlines. The high completeness of the 
galaxy survey enables a detailed investigation of possible en¬ 
vironmental dependence of extended gas around galaxies. To 
increase the number of massive galaxies in the study, we in¬ 
clude eleven galaxies with stellar masses of log M^/Mq >11 
from the Sloan Digital Sky Survey (SPSS; [York et al.|2000 | 
probed in absorption with archival COS quasar spectra. We 


refer to these galaxy samples collectively as the extended- 
CGM (eCGM) survey and combine it with 42 galaxies from 
the COS-Halos program | Werk et al.|2012 Tumlinson et al. 
20131 which are primarily at d < 150 kpc. Together, the two 


samples contain 148 galaxies and enable an investigation of 
the absorption properties of galaxy halos over two decades 
in projected distance. 

The paper proceeds as follows: In Sectionj^we describe 
the galaxy sample. In Section]^ we describe the correspond¬ 
ing absorption-line measurements. In Section]^ we charac¬ 
terize the observed HI and O VI absorption as a function 
of projected distance. In Section we discuss the implica¬ 
tions of our survey results. Throughout the paper, we adopt 
a A cosmology with Dm = 0.3, Da = 0.7, and Hq = 70 
kms“^ Mpc“^. 


2 THE EXTENDED-CGM GALAXY SAMPLE 

To assemble the eCGM galaxy sample, we combined 
our own absorption-blind survey data targeting galaxies of 
Tab < 23 mag in the helds of HE 0226-4110, PKS0405- 
123, LBQS 1435-0134, and PC 1522-1-101 with spectroscopic 
galaxies in the Sloan Digital Sky Survey Data Release 10 
(SDSS; |York et al.||2000| |Ahn et al.||2014|). The HE0226- 
4110, PKS 0405-123, LBQS 1435-0134, and PG1522-kl01 
helds were selected because of the high completeness lev¬ 
els (> 80%) achieved by our surveys for galaxies as faint as 
L = 0.1 L* at a < 0.4 and d < 500 kpc. At smaller projected 
distances of d = 100 and 200 kpc, the survey complete- 
ness increases to 10 0% and 90% respectively (see Figure 2 of 
2013| P| We included SDSS galaxies to increase 


Johnson et al. 


the number of eCGM sample members of log M*/M 0 > 11. 
The resulting galaxy sample is summarized in Table ^ and 
its constrnction is described in this section. 


starburst galaxies ( [Borthakur et al.|2013| |, and more massive lu¬ 
minous red galaxies exhibit reduced Mg II absorption ( [Gauthi^ 
et al.|2010 1. 


^ Here, we adopt a characteristic luminosity of L* galaxies of 
Mr = —21.5 based on the luminosity function measurements from 


Montero-Dorta & Prada (2009 i and Loveday et al. i 2012 I. 























































3 



-16 -18 -20 -22 -24 

[mag] 

Figure 1. Residuals between rest-frame g- and r-band based stel¬ 
lar masses estimated using Equations]^ andand those from the 
NASA-Sloan atlas versus absolute r-band magnitude, Mr- Each 
point represents a NASA-Sloan galaxy. The dashed red line shows 
the median residuals in Mr bins with widths of 0.25 mags. The 
dotted red lines mark zblcr from the median in each bin. Within 
the NASA-Sloan atlas, the systematic error is less than 0.02 dex 
with la scatter of less than 0.15 dex. 


Our galaxy surveys in the HE 0226-4110, PKS 0405-123, 
LBQS 1435-0134, and PG 1522-1-101 fields were carried out 


with the IMACS (Dressier et al. 2011 

1 and LDSS3 spec- 

trographs on the Magellan Telescopes ( 

Shectman & Johns 


20031 targeting galaxies as far as A0 = 10' from the quasar 


sightline (corresponding to d = 1.1 Mpc at 2 = 0.1). The re¬ 
sulting redshift catalogs and corresponding absorption-line 
measurements for the HE 0226-4110, PKS 0405-123 fields 


were published in Chen & Mulchaey (20091 and Johnson 


et al. (20131. The published survey data include deep B-, 


V-, and i?-band photometry in the HE 0226-4110 field. The 
LBQS 1435-0134, and PG 1522-1-101 galaxy redshift catalogs 
and corresponding IMACS g-, r-, and i-band photometry 
will be published in Johnson et al. in prep. For the PKS 0405- 
123 field, we retrieved g-, r-, and i-band images from the 
Mosaic H CCD Imager on the Blanco 4-m telescope from 
the NOAO archive (PI: Brian Keeney; PID=2008B-0194). 
The Mosaic H images consist of 7 x 75 second exposures 
in each bandpass. We measured the galaxy photometry us¬ 
ing Source Extractor | Bertin fc Arnouts|1996 1 with settings 
described in |Johnson et al.| | |2014p . 

To increase the number of massive galaxies in the eCGM 
sample, we cross-matched galaxies of log M^/Mq > 11 in 
the SDSS DRIO spectroscopic survey with quasar sightlines 
that have public, high signal-to-noise COS spectra published 


in Danforth et al. (20141. For the SDSS galaxies, we adopted 
U-, g-, r-, i-, and 2-band cmodel magnitude^ 


From this parent sample of galaxies, we selected those 
with projected distances from the quasar sightline d < 1 
Mpc, within the redshift range z ~ 0.1 — 0.4, and with red¬ 
shift at least 10, 000 kms“^ below the quasar redshift. O VI 
detections have been reported at projected distances of up to 
300 kpc from sub-L* galaxies ( Prochaska et al.|2011 1 which 
corresponds to 2 — 3 times the galaxy halo virial radius. 
The d = 1 Mpc cut was chosen to be sufficiently large to 
search for absorption at d/Rh = 2 — 3 around more massive 
galaxies that have virial radii of i?h ~ 300 kpc. The red¬ 
shift lower limit was chosen so that the O VI doublet falls 
within the COS wavelength range. We note that the red¬ 
shift range for which the OVI doublet is observable with 
COS depends on the settings used for the observations. The 
redshift lower limit therefore varied slightly from sightline 
to sightline. The redshift upper limit was motivated by the 
depth of our IMACS redshift surveys which are sensitive to 
galaxies of L > 0.1 L* at 2 < 0.4. The requirement that the 
galaxy redshift be at least 10,000 kms“^ below the quasar 
redshift was imposed to avoid confusion with gas outflowing 
from the quasar (e.g. Wild et al.|2008 K In total, 106 galaxies 
met these criteria with eleven galaxies from the SDSS and 
95 from our survey data. 

For each galaxy, we measured the rest-frame, g- and 
r-band absolute magnitudes using the multi-band optical 
photometr^p] described above and the kcorrect tool (Blan 


ton fc Roweis|2007 We then estimated mass-to-light ratios 


using the relation shown in Equation which is a fit of the 
stellar mass to r-band light ratio (in solar units) for low 


redshift galaxies in the NASA-Sloan atlas (e.g. Mailer et al. 
|2009p . For galaxies redder than Mg — Mr = 0.15 in the 
rest-frame, the fit is a quadratic in Mg — Mr- For galaxies 
bluer than Mg — Mr = 0.15, we set the mass-to-light ratio 
to log (M,/I/r) = —0.6 which corresponds to an imposed 
minimum age of « 0.3 Gyr for a simple stellar population 
with solar metallicity. 


i -0.6 Mg-Mr<0.lh 

-1.01 + 2.95(Mg-MQ... Mg-M,.>0.15 (1) 

-1.67(Mg - Mrf 

We then estimated the stellar mass of the eCGM galaxies 
using Equationj^and the mass-to-light ratios from Equation 

m 

log {M,/MQ)gr = 1.872 - 0.4M,. -b log {M^/U) (2) 


The rest-frame g- and r- band stellar mass estimates 
from Equations and reproduce the NASA-Sloan stel¬ 
lar masses which are calculated with kcorrect, a |Chabrier| 
120031 initial mass function and GALEX FUV + NUV and 
SDSS ugriz photometry with a systematic error of less than 
0.02 dex and la scatter of less than 0.15 dex for galaxies of 
Mr = —16 to —24 as shown in Figure The NASA-Sloan 
atlas contains low-redshift galaxies of 2 < 0.055 so the use 
of a mass-to-light ratio relation based on these galaxies for 
the 2 = 0.1 — 0.4 eCGM sample may result in systematic 
errors due to galaxy evolution over cosmic time. To esti¬ 
mate the magnitude of this systematic error, we compared 
stellar masses calculated using Equations and to those 


® The results presented in this paper do not change if model 
rather than cmodel magnitudes are adopted for SDSS galaxies. 


Corrected for Milky Way extinction from Schlegel et al. 


119981. 
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Figure 2. Distribution of projected distance, d, between each 
galaxy and the nearby COS quasar sightline versus stellar mass 
for the eCGM galaxy sample (black circles) and the COS-Halos 
survey (magenta squares). The top three panels show the stellar 
mass, projected distance, and redshift distributions of the eCGM 
and COS-Halos galaxy samples. 


based on precise ugriz photometry for 2 = 0.1 — 0.4 galaxies 
from the PRIMUS galaxy survey ( Coil et al.|2011 1. The sys¬ 
tematic error introduced by using a 2 < 0.055 relation for 
2 = 0.1 —0.4 galaxies is less than 0.1 dex over the luminosity 
range spanned by the eCGM sample. 

The eCGM galaxy sample spans a stellar mass range of 
log M^/Mq = 8.4 — 11.5 with a median of log M^/Mq = 
10.3 as shown in Figure 

In order to compare absorption properties of galaxies 
across a range of stellar mass, we measured projected dis¬ 
tance in units of the galaxy dark matter halo virial radius 


defined by Bryan & Norman ( 1998[) using the stellar-to-halo 


mass relation from Kravtsov et al. (20141. To ensure consis¬ 


tency across the sample, we re-measured the virial radii of 
the COS-Halos galaxies using the same methods. We note 
that our virial radius estimates are as much as 50% smaller 


than those of Tumlinson et al. (20131 because of the use 


of a different stellar-to-halo mass relation. The differences 
between the stellar-to-hallo mass relations from IKravtso^ 


et al. (20141 and previous work such as Moster et al. (20101 


and Behroozi et al. (20131 are driven by systematic errors 


in the galaxy photometry in the catalogs used by |Moster 


eCGM galaxy sample to isolated galaxies only. Instead, we 
classified galaxies as “non-isolated” or “isolated” based on 
the presence or absence of nearby neighbors. In the HE 0226- 
4110, PKS 0405-123, LBQS 1435-0134, and PC 1522-blOl 
fields, where our galaxy survey data are highly complete, 
we classified galaxies as “non-isolated” if they have a spec¬ 
troscopic neighbor within a projected distance of less than 
500 kpc, a radial velocity difference of |Au| < 300 kms“^, 
and with stellar mass of at least one-third of that of the 
eCGG or COS-Halos survey member. Otherwise, we classi¬ 
fied galaxies without such nearby neighbors as “isolated”. 

Among non-isolated galaxies, the relevant projected dis¬ 
tance between the quasar sightline and the galaxy is not 
clear apriori, and care must be taken in defining the distance 
to a non-isolated galaxy to prevent a bias when comparing 
with isolated ones. For example, a mass or light-weighted 
projected distance could be large even though the gas de¬ 
tected in absorption may be more closely associated with 
a less massive “group” member close to the quasar sight¬ 
line. To avoid a bias that could make the gas around non¬ 
isolated galaxies appear artificially extended, we define the 
distance to non-isolated galaxies as the minimum value of 
d/Rh among “group” members. This distance definition is 
a self-consistent choice for comparing the spatial extent of 
absorbing gas around isolated and non-isolated galaxies. 

The available galaxy survey data in the COS-Halos and 
SDSS fields are not sufficiently complete to confidently clas¬ 
sify galaxies by environment based on spectroscopy alone. 
To classify these galaxies, we performed a literature search 
to identify members of previously known galaxy groups (e.g. 


McConnachie et al. 20091 and clusters (e.g. Koester et al. 


20071. In addition, we visually inspected the SDSS images 


in the vicinity of each galaxy to search for likely neighbors 
based on galaxy brightness, size, color, and SDSS photo¬ 
metric redshift. We classified galaxies with such neighbors 
as non-isolated and those without as isolated. 


3 ABSORPTION DATA 

For each galaxy in the eCGM sample, we searched for 
HI and O VI absorption in the COS quasar spectra within a 
radial velocity interval of Av = ±300 kms“^ of the galaxy 
systemic redshift. Tumlinson et al. (20111 found that the 


radial velocity distribution of O VI absorption-components 
around galaxies at d < 150 kpc is characterized by a stan¬ 
dard deviation of cr « 100 kms“^. The Av = ±300 kms“^ 
search window was chosen to include the vast majority of ab¬ 
sorption associated with the eCGM and COS-Halos galax- 
ie^ If HI or OVI absorption was detected, we adopted 
available Voigt profile fitting results from [Johnson et al. 


(|2013[), Tumlinson et al. (20131, Werk et al. (20131, and 


Sav- 


et al. (20101 and Behroozi et al. (20131 (see 

Bernardi et al. 

age et al. (2014). When these measurements are not avail- 

2013). The stellar-to-halo mass relation from 

Kravtsov et al. 

able, we performed our own Voigt profile fitting using the 


(20141 is in good agreement with independent constraints 
from weak leasing (e.g. Reyes et al.|2012 l. 

To search for correlations between CGM properties and 
the presence of star-format ion, we spectroscopically clas¬ 
sified galaxies as late-type (emission-line dominated) and 
early-type (absorption-line dominated) based on a PCA fit¬ 
ting technique described in Chen & Mulchaey ( |2009| . 

Unlike some previous surveys, we did not restrict the 


son et al. (20131. If HI or O VI absorption was not detected, 
we placed 2a upper limits on the equivalent width, inte¬ 
grating over a velocity interval of 75 kms“^ for HI and 90 


® The results presented in this paper do not change if the velocity 
window is increased to Av = ±600 kms“^. 
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d/R^ d/R^ d/R^ 


Figure 3. Column density within An = ±300 kms“^ versus projected distance measured in units of the galaxy virial radius with HI 
in the top panels and O VI on the bottom. The full galaxy samples are shown in the left column and are separated into isolated and 
non-isolated members in the middle and right columns respectively. Late- and early-type galaxies are shown in blue and red respectively. 
eCGM galaxies are displayed as filled circles and COS-Halos galaxies as open squares. Non-detections are shown as 2ct upper limits 
marked by downward arrows and lighter coloring. The vertical error bars represent Icr statistical uncertainties reported by VPFIT and 
do not include uncertainty due to continuum placement. Galaxies with upper limits greater than logV/cm"^ = 14.0 are not considered 
to have sensitive constraints on the g as column density and are not displayed. The mean absorption expected from random sightlines 


based on the 


d?j\r 


measurement from 


Danforth et al. 


(20141 is shown as horizontal dotted line. 


kms“^ for OVI. We then converted this equivalent width 
limit to a column density assuming that the gas is opti¬ 
cally thin. The 75 kms“^ and 90 kms“^ integration win¬ 
dows correspond to the median full-width-at-half-maximum 
of unsaturated HI and OVI absorption components found 
in the eCGM survey. 


For 19 galaxies, [Tumlinson et ah] ( |2013[ ) report only 
lower limits on HI column density due to saturation in all 
available Lyman series transitions. For these galaxies, we 
place upper limits on the possible HI column density based 
on the lack of damping wings detected in Lya. To do so, 
we measured the Lya absorption equivalent width and con¬ 
verted this to a HI column density upper limit assuming a 
single component with doppler width 6 = 20 kms“^. 


4 THE EXTENT OF HI AND O VI 

ABSORPTION AROUND GALAXIES 


With the galaxy and absorption data described in Sec¬ 
tions and we characterize the HI and O VI column 
densities of the eCGM and COS-Halos galaxy sightlines 
as a function of projected distance in Figure The cov¬ 
ering fractions for strong HI and OVI absorption sys¬ 
tems of logA(HI)/cm-2 > 14.0 and log A(0 Vlj/cm’^ > 
13.5 are shown in Figure The uncertainties in cover¬ 
ing fractions are calculated using standard binomial statis¬ 
tics, and we do not include galaxies with non-detections 
with upper limits greater than log Al(HI)/cm~^ > 13.5 or 
log A(0 VI)/cm“^ > 14.0 in the calculations. When multi¬ 
ple galaxies are associated with the same absorption system, 
we show only the galaxy at the smallest d/Rh. 
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d/R^ d/R^ d/R^ 


Figure 4. Covering fraction for strong HI and O VI absorption systems as a function of d/Ri^. The top panels show the covering fraction 
for HI column density greater than log _/Vo(HI)/cm“^ = 14.0 and the bottom panels show the covering fraction for O VI column density 
greater than log Alo(0 VI)/cm“^ = 13.5. The late-type galaxy samples are shown in blue and early-type in red. The full eCGM and 
COS-Halos samples are shown on the left, isolated galaxies in the middle, and non-isolated ones on the right. The horizontal error bars 
mark the full range of galaxies within each bin and vertical error bars represent 68% confidence intervals calculated using standard 
binomial statistics. We do not include non-detections with column density limits greater than Nq in the covering frac tion calculation. 
For both HI and O VI, we show the expected covering fraction for random sightlines based 
| |2014| |. 


measurements from 


Danforth et al. 


4.1 HI 


Consistent with previous surveys (e.g. Chen et al.|2001[ |, 
we find a clear trend of decreasing HI column density with 
increasing projected distance (see the top panels of Figures]^ 
and|^. Lyman-limit systems (LLS) and partial Lyman-limit 
systems (log A'^(HI)/cm“^ > 16) are common at d < i?h, 
but none are found at larger distances. Considering the 
eCGM and COS-Halos samples together, we identify strong 
HI systems with log V(HI)/cm“^ > 14 around 48 of 54 
galaxies at d < i?h indicating a mean covering fraction of 
(«Hi) = 0.89to'o6' At larger distances of d/Rh = 1 — 3 such 
HI absorption systems are found for only 7 of 43 galax¬ 
ies ((khi) = O.IGIq oI). The HI covering fraction decreases 
to levels consistent with expectations from coincidental ab¬ 
sorption systems unrelated to the galaxy survey members at 
d > 3 Rh ■ 

To determine whether galaxy environment has an influ¬ 


ence on HI absorption far from galaxies, we show the column 
densities and covering fractions for isolated and non-isolated 
galaxies in the bottom middle and bottom right panels of 
Figures 1^ and 1^ Though the only two galaxies detected at 
d > with log V^(HI)/cm“^ > 15 have nearby neigh¬ 
bors, there is no clear enhancement in HI absorption in the 
non-isolated sample. A series of logrank tests over the range 
of d/i?h = 1 — 3 finds no evidence for differential HI ab¬ 
sorption between the isolated and non-isolated samples. We 
note, however, that contamination from coincidental Lya 
absorbers in the IGM is non-neglegible. This contamination 
may obscure a difference in outer-halo HI absorption be¬ 
tween the two samples. 

Cross-correlation studies such as |Chen fc Mulchaey| 


(20091 and Tejos et al. (2014l indicate that strong HI ab¬ 


sorption systems primarily reside in the gaseous halos of 
star-forming galaxies. To investigate the possibility of en¬ 
hanced HI absorption in galaxies with recent star-formation 
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activity, we show late- and early-type galaxies in blue and 
red respectively in Figures and Late-type galaxies ex¬ 
hibit excess HI absorption relative to early-type galaxies at 
projected distances of d < _Rh- In particular, we find that 
35 of 36 late-type galaxies at d < Rh are associated with 
HI systems of log A^(HI)/cm~^ > 14.0 indicating a mean 
covering fraction of (khi) = 0.97lg gg compared to 13 of 
18 early-type galaxies ((khi) = 0.72to;?2)- However, a com¬ 
parison of the distributions of HI column densities found 
at d < i?h around late- and early-type galaxies is compli¬ 
cated by the poorly determined HI column densities for the 
LLS and partial-LLS that are common at small projected 
distances (see discussion in Tumlinson et al.]|2013 |. 

At larger distances of d/Rh = 1 — 3 the HI covering 
fraction for systems with log Af(HI)/cm“^ > 14 are simi¬ 
lar for late- and early-type galaxies. However, covering frac¬ 
tions are somewhat insensitive to the underlying HI col¬ 
umn density distribution. To compare the HI column den¬ 
sity distribution of late- and early-type galaxies in the outer 
halo, we perform a logrank test and find that the probability 
that the column densities around the two galaxy classes at 
d/i?h = 1 — 3 are drawn from the same underlying distribu¬ 
tion to be P = 2%. To investigate the nature of this likely 
differential HI absorption around late- and early-type galax¬ 
ies at large projected distances, we compare the HI column 
densities of the two samples using the Kaplan-Meier estima¬ 
tor, an unbiased, non-parametric estimator of cumulative 


distributions in the presence of upper limits (see Feigel- 
son fc Nelson|19 ^ \ in Figure]^ We find excess HI absorp¬ 


tion around late-type galaxies driven primarily by moderate 
strength absorption systems of log ^^(Hlj/cm”^ « 13.5 to 
14.0. Based on the Kaplan-Meier curves, we infer a median 
HI column density of (log A'(HI)/cm“^) rs 13.8 for late- 
type galaxies at d/Ph = 1 — 3. For early-type galaxies, the 
fraction of HI detections at these distances is too small for 
a robust measurement of the median, but we are able to 
place an upper limit on the median HI column density of 
(log iV(HI)/cm-2) < 13.4. 



Figure 5. Cumulative fraction, P, of galaxies at d/R\^ = 1 — 3 
with total HI column densities no greater than log Al(HI)/cm“^. 
The late-type galaxies are shown in blue solid line and early-type 
as red dashed line. The cumulative distributions are estimated 
using the Kaplan-Meier curve, an unbiased, non-parametric es¬ 
timator of the cumulative distribution in the presence of upper- 
limits. The Kaplan-Meier estimator changes value at detections 
but not at non-detections. For clarity, non-detections with upper 
limits are shown as open symbols. The shaded bands represent 
68% confidence intervals including uncertainty due to both sam¬ 
ple variance and column density measurement errors calculated 
with combined boot-strap and Monte-Carlo resampling. At col¬ 
umn densities smaller than the highest upper limit that is less 
than the lowest detection, the Kaplan-Meier estimator represents 
an upper limit on the cumulative distribution. These column den¬ 
sity ranges are shown in faded color. 


4.2 OVI 

Like HI, we find a clear anti-correlation between 
OVI column density and projected distance (see the bot¬ 
tom left panel of Figure [^. Strong OVI systems with 
log V(0 VI)/cm“^ « 14.5 — 15.0 are common at d < i?h 
but are rare at larger projected distances. O VI systems 
with log ^"(0 VI)/cm“^ « 13.5 — 14.0 are detected out 
to d/i?h = 2.3, but at larger distances, no such strong 
OVI systems are found. Considering the eCGM and COS- 
Halos samples as a whole, we find OVI absorption with 
log V(0 VI)/cm“^ > 13.5 for 36 of 42 galaxies at d < i?h 
indicating a mean covering fraction of (kqvi) = 0 . 86 lQ g 7 . 
Isolated galaxies exhibit OVI absorption out to d ~ f?h, be¬ 
yond the “metal-boundary” seen in CIV and lower ions at 
d « 0.7 dZh- At d/iZh = 1 — 3 strong OVI absorption systems 
are found for only 5 of 34 galaxies ((kovi) = 0.157Q g|). 
At larger distances of d > 3iZh, the OVI covering fraction 
decreases to levels consistent with expectations from coin¬ 
cidental absorption systems unrelated to the galaxy survey 
members (see the bottom left panel of Figure]^. 

The high completeness of the galaxy survey enables a 
detailed investigation of possible environmental dependence 


of extended gas around galaxies, and we show O VI column 
density and covering fraction versus d/iZh for isolated and 
non-isolated galaxies separately in the bottom middle and 
bottom right panels of Figures and The decrease in 
incidence of strong OVI systems at d « iZh is particularly 
sharp among isolated galaxies. Of the 51 isolated galaxies 
at d > iZh in the sample we detect log A(0 VI)/cm~^ > 
13.5 absorption systems for only two galaxies (at d/7Zh = 
1.1 and 1 . 2 ) indicating a mean covering fraction of kqvi = 
0.04(ig g7 In contrast, three non-isolated galaxies at d/iZh = 
1.3,1.8, and 2.3 are found to have absorption systems with 
log N (O VI)/cm“^ > 13.5 out of a sample of 24 such galaxies 
at d > iZh (kovi — 0.13lo (, 4 ). This suggests that the OVI 
absorption around non-isolated galaxies is more extended 
than around isolated galaxies. 

To evaluate the possibility that non-isolated galaxies 
have more extended O VI absorbing gas, we compare the 
OVI column density distributions of the isolated and non¬ 
isolated galaxy samples at d > 1.5 iZh. A logrank test finds 
that the probability that the O VI column densities of the 
isolated and non-isolated galaxies at d > 1.5 iZh are drawn 
from the same underlying distribution to be P < 1%. The 
results of the logrank test are insensitive to the choice of 














d > 1.5-Rh cut-off and return a result of P < 1% over the 
full range of cut-off distances between 1.2 and 1.8 Ph and 
P < 2% between 1.1 and 2.2Ph- These findings indicate 
that galaxies with nearby neighbors exhibit more extended 
O VI absorbing gas than isolated galaxies. 

Though the completeness of our galaxy survey is high, 
it is possible that the detections of O VI at large distances 
from non-isolated galaxies could be due to the presence of 
an additional neighbor closer to the sightline that was not 
observed in our survey. However, the survey incompleteness 
applies to the isolated and non-isolated galaxies equally. In 
addition, the galaxy survey completeness levels increase at 
small projected distances from the quasar sightline (see Fig¬ 
ure 2 of Johnson et al.|2013 |. In particular, at the redshifts 
and distances of the non-isolated galaxies detected in O VI 
absorption at d > Ph, our galaxy survey data are > 90% 
complete for galaxies as faint as 0.1 P*. 

Previous surveys (e.g. |Chen fc Mulchaey||2009l |Tum-| 
linson et al. 20111 found that late-type galaxies exhibit 


enhanced OVI absorption relative to more massive early- 
type galaxies at d < 150 kpc. We confirm that late-type 
galaxies show enhanced OVI absorption at d < Ph rel¬ 
ative to early-type ones. However, OVI absorption sys¬ 
tems with log V(0 VI)/cm“^ > 13.5 are found for six of 
eleven early-type galaxies at d < Ph indicating a non- 
negligible covering fraction of (kovi) ~ 0 .55lo'i4. At larger 
distances of d/Ph = 1 — 3 we find OVI absorption with 
log V(0 VI)/cm“^ > 13.5 for four of 14 late-type galax¬ 
ies (kovi) ~ O. 29 I 0 Jg compared to one of 19 early-type 
galaxies (kovi) ~ O.OdlQ'J®. Interestingly, the only early- 
type galaxy exhibiting a strong OVI system at d > Ph is 


the brightest member of a maxBCG cluster (Koester et al. 


2007) at d = 1.3 Ph from the quasar sightline. A logrank 


test finds that the probability that the OVI column den¬ 
sities found at d/Ph = 1 — 3 from late- and early-type 
galaxies are drawn from the same underlying distribution 
to be P = 1%. We note that the O VI excess in the outer 
halos of late-type galaxies is driven by absorption systems 
with log V(0 VI)/cm“^ « 13.5 — 14.0 while the excess at 
smaller projected distances is driven by stronger systems 
with log V(0 VI)/cm“^ « 14.5. 


5.1 Possible environmental effect in distribnting 
heavy elements to d > Ph 

Among isolated galaxies without neighbors within 500 
kpc and with stellar mass ratio greater than one-third, we 
find that OVI absorption is confined to within d < 1.2Ph. 
None of the 18 isolated galaxies at d/Ph = 1.2 — 3.0 in 
our survey are found to have O VI absorption in the COS 
quasar spectrum. In contrast, three out of ten galaxies with 
nearby neighbors probed at d/Ph = 1.2 — 3.0 are found to 
exhibit OVI absorption with log V(0 VI)/cm“^ > 13.5. On 
the other hand, no evidence for differential HI absorption is 
found between the isolated and non-isolated galaxy samples. 
These findings suggest that galaxy interactions play a key 
role in distributing heavy elements to large distances from 
galaxy centers, well beyond the enriched gaseous halos that 
surround individual galaxies. 

Galaxy interactions can produce heavy elements at 
large projected distances through tidal stripping during 


satellite accretion (e.g. Wang 19931 or ram-pressure strip¬ 
ping (e.g. Gunn & Gott 1972). Recently, HST images of 
galaxy clusters have revealed that ram pressure stripping 
can remove and compress copious amounts of gas from in¬ 
falling galaxies as evidenced by vigorous starbursts in debris 
trails ( Ebeling et al.|2014 ). Moreover, a galaxy survey in the 
field of an ultra-strong Mg H absorption system near a lumi¬ 
nous red galaxy revealed a group of evolved galaxies with no 
evidence of recent star-formation activity ( Gauthi^|2013 ). 
The lack of star-formation in the group and kinematics of 
the Mg II absorption are most readily explained if the Mg H 
absorber arises from stripped gas. Finally, in [Johnson et ah] 
| |2014[ ) we presented the discovery of a “transparent” sight¬ 
line with no strong H I, Mg H or O VI absorption at d < 20 
kpc from a pair of strongly interacting galaxies separated 
by a projected distance of 9 kpc. The lack of strong absorp¬ 
tion systems detected at such small projected distance from 
this galaxy pair can be explained if the inner halo gas of the 
two galaxies has been stripped to larger distances during the 
galaxy interaction. 


5 DISCUSSION 

Using a highly complete survey (> 80%) of faint galax¬ 
ies (L > 0.11/*) at 2 < 0.4 and d < 500 kpc in the fields of 
four GOS quasar sightlines with high signal-to-noise spectra 
available in the HST archive, we have searched for the key 
features that determine the extended gaseous properties of 
galaxy halos in HI and O VI absorption. Our main findings 
are: 

(i) galaxies with nearby neighbors exhibit enhanced O VI 
absorption at d > Ph relative to isolated galaxies; 

(ii) among isolated galaxies, O VI absorption extends to 
d ~ Ph, beyond the extent observed in lower ionization 
species such as Si III and GIV; and 

(iii) late-type galaxies exhibit enhanced HI and O VI ab¬ 
sorption beyond characteristic halo radii at d/Ph = 1 — 3. 


5.2 Comparison with Prochaska et al. 2011 


Prochaska et al. (20111 found that galaxies at 2 ; < 0.1 


and of —21.1 < Mr < —18.6 exhibit near unity covering 
fraction at d < 300 kpc and concluded that such “sub-L*” 
galaxies possess OVI bearing gaseous halos that extend to 
d « 3Ph. Our highly complete galaxy survey data shed 
new light on these observations, indicating that O VI ab¬ 
sorption at such large distances from galaxy halos arise in 
systems of multiple galaxies. Indeed, a literature search of 
the NASA Extragalactic database for spectroscopic neigh¬ 
bors reveals that all of the galaxies of —21.1 < Mr < —18.6 
from |Prochaska et al.| ( [20TT| with O VI absorption at d > Ph 
have neighbors within d = 500 kpc and satisfy our defini¬ 
tion of a non-isolated galaxy (see Figure]^. The data from 
[Prochaska et ah] ( [2011[ ) are therefore consistent with our 
conclusion that isolated galaxies exhibit low covering frac¬ 
tions at d > 1.2 Ph and support the possibility that galaxy 
interactions are effective at producing O VI absorbing gas at 
large projected distances. 
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Figure 6. O VI column density versus dj Hh for sub-L* galaxies of 
—21.1 < Mr < —18.6. Galaxies from the eCGM survey are shown 
as black circles, those from|Werk et al.|l[2012| l as magenta squares, 
and those from |Prochaska et al.| ( |2011^ as orange triangles. Non¬ 
isolated galaxies are outlined by large black symbols. 


5.3 OVI absorption in isolated galaxies beyond 
the “metal-boundary” at d/i?h = 0.7 

Even among isolated galaxies, O VI absorption extends 
beyond the “metal-boundary” observed by [Liang fc Chen| 
(20141 at d = 0.7 i?h in lower-ionization state, enriched 


gas traced by absorbers such as Sill, Si III, CII, and CIV. 
Specifically, we find OVI absorption for seven of ten iso¬ 
lated galaxies at d/Rh = 0.7 to 1.2. The presence of OVI 
absorption in the outer regions of isolated galaxy halos indi¬ 
cates some heavy element enrichment at large radii and that 
highly ionized, OVI bearing clouds can form or survive in 
the outer halo while Si II, Si III, and CIV bearing clouds can¬ 
not. Using the mean absorption found in a stack of isolated 
galaxies in the eCGM and COS-Halos samples at d/Rh = 
0.7 — 1.2, we estimate the mean metallicity of these clouds. 
In the stack, we measure absorption of log V(H I)/cm“^ = 
14.4, log V(0 VI)/cm“^ = 13.9, and place an upper limit 
on cm absorption of log V(CIII)/cm“^ < 12.6. Assuming 
photoionization equilibrium models calculated with Cloudy 
(Ferland et al. 1998[) Version 1 0.0, a photoionization back- 
ground from jHaardt &: Madauj ( |2012[ ) at z = 0.2, and solar 
abundance patterns, these mean column densities are consis¬ 
tent with gas in photoionization equilibrium with ionization 
parameter log U = —1.1 and a metallicity of one-tenth so¬ 
lar. The OVI and HI bearing gas may not be co-spatial 
allowing the ionization parameter of the HI gas to be lower. 
If this is the case, then the upper limit on C III absorption 
places stricter upper limits on the mean gas metallicity. We 
therefore place an upper limit of one-tenth solar on the mean 
metallicity of outer-halo HI gas clouds at d « i7h. This is 
consistent with the metallicity upper limit of one-tenth solar 


found by Liang & Chen (2014 1 based on CIV non-detections 


at d > 0.8 i7h. 


5.4 Differential HI and OVI absorption between 
late- versus early-type galaxies 

We find excess HI and O VI absorption around late- 
type galaxies relative to early-type ones at d/i?h = 1 — 3 with 
99% and 98% confidence respectively, extending previous re¬ 


sults from Chen & Mulchaey (20091, Tumlinson et al. (2011 1 , 


and Tejos et al. (20141 to larger distances. The OVI excess 


around late-type galaxies could be explained if the ioniza¬ 
tion mechanisms that produce O VI absorbing gas are the 
result of star-formation (e.g. photoionization from a young 
stellar population). However, this is not consistent with the 
observed HI excess unless the ionization mechanism is effec¬ 
tive at producing highly ionized heavy element species but 
ineffective at ionizing hydrogen (e.g. Auger ionization). 

Higher levels of HI and O VI absorption around late- 
type galaxies could be simultaneously explained if the ab¬ 
sorption systems trace starburst driven outflows. However, 


simulations from Oppenheimer & Dave (2008 1 that incor¬ 


porate strong, momentum-driven winds find that these out¬ 
flows are characterized by a turnaround radius of i?turn = 
80 ±20 kpc with only weak dependence on mass and redshift 
at z < 1. For the late-type galaxies in our survey, this corre¬ 
sponds to only i7turn/7?h ~ 0.4, far less than ci/i?h = 1 — 3. 
Alternatively, the excess HI and O VI absorption at large 
distances from late-type galaxies can be explained if HI, 
OVI, and star-formation are observational signatures of ha¬ 
los in cosmic environments with multi-phase gaseous reser¬ 
voirs. 

The late-type members of the eCGM and COS- 
Halos galaxy samples are systematically less massive than 
the early-type members with mean stellar masses of 
log(M*/M 0 ) = 10.2,10.9 corresponding to inferred halo 
masses of log {Mh/Mq) = 11.7,12.6 respectively. The mean 
doppler parameter of the HI absorption components de¬ 
tected at d/7?h = 1 — 3 from eCGM galaxies is fo = 40 
kms“^ which indicates that the HI absorbers trace cool- 
warm, T < 10® K gas. Our survey data are therefore con¬ 
sistent with a decreased incidence of cool-warm gas around 
more massive halos. This extends previous results from jGau^ 
jthier et al.| ( | 2010 | ) and |Yoon et al.| ( |2012[ ) who found reduced 
cool-warm gas absorption in > 10 uW© mass halos relative 
to less massive ones. A larger sample of early-type galaxies 
probed in absorption is required to differentiate the effects 
of recent star-formation activity and mass. 
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